Giant mass and anomalous mobility of particles in a fermionic system 
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We calculate the mobility of a heavy particle coupled to a Fermi sea within a non-perturbative 
approach valid at all temperatures. The interplay of particle recoil and of strong coupling effects, 
leading to the orthogonality catastrophe for an infinitely heavy particle, is carefully taken into 
account. We find two novel types of strong coupling effects: a new low energy scale T* and a giant 
mass renormalization in the case of either near-resonant scattering or a large transport cross section 
a. The mobility is shown to obey two different power laws below and above T* . For u > Aj, where 
Xf is the Fermi wave length, an exponentially large effective mass suppresses the mobility. 



Some remarkable advances in condensed matter 
physics are based on a thorough analysis of models of a 
single impurity coupled to a many-particle system. The 
polaron problem and the mobility of ions in liquid 3 He 
are well-known examples for the dynamical behavior of 
a delocalized particle in a bosonic or fermionic bath, re- 
spectively. We will consider the latter case of a delocal- 
ized but heavy particle strongly interacting with a Fermi 
liquid. 

Concerning the experimental context, a 3 He + -ion in 
normal-fluid 3 He is the cleanest possible heavy-particle 
Fermi-liquid system. Unfortunately, the transition to the 
superfluid state cuts off the available low-temperature 
range jj] . The particle is indeed heavy and of large cross 
section since the ion forms a snowball: about 100 He- 
lium atoms (depending on pressure) are tightly bound at 
temperatures T up to several hundred mK. Other real 
fermionic systems, such as muons in metals or valence 
band holes in rt-type doped semiconductors, show a more 
complex dynamical behavior, usually involving several 
diffusion mechanisms Nevertheless, our approach 

is supposed to be valid in the low-T regime. 

We will discuss two strong coupling effects of a heavy 
particle |^J^] in a fermionic bath, both of which have 
not been reported before, basically because a microscopic 
strong coupling calculation of the mobility has not been 
developed for low temperatures. As to the origin of 
strong coupling effects we actually have to distinguish 
two scenarios: (1) the scattering phase at the Fermi en- 
ergy of bath particles is close to S = ir/2 in any single 
scattering channel; (2) the transport cross section a of 
the heavy particle is large: a s> (27r/fc/) 2 . Obviously, 
each of the two scenarios implies a strong coupling of 
bath particles to the impurity. 

In the first scenario we work close to a resonance at 
the Fermi energy. This may signify that an attractive 
channel exists with sufficiently large binding energy. Al- 
ternatively, we may consider a strong enough repulsive 
interaction between the particle and the fermionic bath, 
e.g. the Hubbard model in the spin-polarized sector with 



one overturned spin which interacts with the free Fermi 
system of up-spins through a local interaction U. If the 
tight-binding band is nearly half-filled with up-spins and 
the repulsive interaction approaches infinity, U — > oo, the 
bath particles at the Fermi energy encounter a phase shift 
close to 7r/2. These problems cannot be solved within 
perturbation theory. Indeed, the heavy-particle mobility 
changes its temperature dependence drastically with re- 
spect to the pcrturbational limit if the phase shift exceeds 
a minimal value in at least one channel: the temperature 
(and frequency) power law behavior acquires an anoma- 
lous exponent as we will demonstrate. 

The second scenario may be encountered already for 
relatively small phase shifts in each of the scatter- 
ing channels provided the particle couples to the bath 
through a sufficiently large number N of channels. We 
find that the effective mass of the particle is tremendously 
increased for low T: it scales with an exponential func- 
tion, the exponent of which is proportional to a power of 
N. 

The mobility of a heavy particle in a fermionic environ- 
ment has been investigated for more than three decades 
fjJIIJlll- Concerning translational invariant systems, as 
ions in^He, the work by Wolfle et al. summarized the pre- 
vious findings and presented the first unified treatment 
for all temperature ranges below the Fermi temperature 
within a self-consistent Mori-Zwanzig scheme. The ap- 
proach is phenomenological and cannot account for part 
of the strong coupling effects discussed here. More re- 
cently Prokof'ev || put forward a microscopic path in- 
tegral scheme, yet his evaluation is restricted to the high 
temperature corrections in the mobility. 

The basis of our evaluation (for space dimension d > 1) 
is an effective low-energy action for the heavy particle, 
first introduced by Sols and Guinea investigated in 
more detail by Prokof'ev || , and evaluated for the heavy- 
particle spectral function in a previous paper [jllj : 

S eS = S a -^ J J^dTdT 1 T{\R{t')-R{t)\) P (t'-t). (1) 
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The first term on the rhs is the action for a free par- 
ticle, S = | dr M q R(t) 2 , and T{\R' - R\) is de- 
termined through the overlap of the electronic ground 
state 4>R with particle at position R and a state 4> R > in 
the thermodynamic limit N e —> oo fl^] : (4>r'\4>r) ~~ > 
exp{-[(<5/V) 2 ~T{\R' - R\)} log A^} where 

T(\R\) = (Sin) 2 - ^rarcsin 2 (Vl - x 2 sin 8) (2) 

with x= {e tkR )\k\ = k f (see Fig. |l|). These relations apply 
for s-wave scattering with phase shift 8. Generalizations 
for several scattering channels are found in the literature 
J[2,[L3| . The leading low-energy form of p{r— r') in d > 



I~lf is P (t-t>) = (tt//3)E^oI 



which 



reflects that the number of virtual excitations close to 
the Fermi surface is proportional to uo. This behavior 
and the Friedel oscillations of J-(R) are rooted in the 
fermionic character of the bath. 

Alternatively, one could set up a Hamiltonian for both, 
the heavy particle (P,R) and the fermionic degrees of 
freedom (cj., c*,) 
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for s-wave scattering. In a second step, the fermions 
are traced out. This procedure results in an effective 
heavy-particle action which involves any number of re- 
tardations (i.e. time integrals) and it cannot be handled 
further except for a few limiting cases. Fortunately, the 
effective action (|l|) reproduces all exactly known limits 
of this more basic approach: the spectral function of the 
heavy particle was evaluated and discussed in this re- 
spect and the mobility /J,(T) is identical in both ap- 
proaches in the exactly known high-temperature limit: a 
constant behavior /i(T) = /i Q = e/(M T) plus logarithmic 
corrections P,pT|. Here r = nkja/M is the scattering 
rate, kf — 2n/Xf the Fermi wave vector and n the parti- 
cle density. 

The success of the approach based on the effective 
action ([!]) is actually not surprising. It correctly ac- 
counts for the two competing effects which control the 
dynamics: the orthogonality catastrophe (OC) through 
the exact form of the overlap and the recoil through 
translational invariance. The relevant low-energy scale 
is supposed to settle which of the two effects dominates. 
It may be anticipated that the maximum recoil energy 
Er = (2kf) 2 /2M Q transferred in a collision process is the 
relevant energy scale. For frequency to or temperature T 
far above Er the energy balance in scattering processes 
is independent of the recoil — i.e. physics is governed by 
the R = saddle point and the corresponding fluctua- 
tions: the high-T mobility results from the coefficient of 
the i? 2 -term in T{R) = (8/tt) 2 - (TM /2ir) R 2 + 0(R 4 ). 
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FIG. 1. Qualitative behavior of T(R) for strong coupling. 

For frequencies or temperatures much less than Er 
the available phase space for scattering processes with 
momentum transfer |Afe| ~ 2k / is restricted by recoil, 
i.e. in d dimensions by a factor (T/Er)^ -1 ^ 2 and we 
expect a well-defined quasi particle for T — > — albeit 
with renormalized mass or possibly coupling constant. 
For the mobility this corresponds to a 1/T^ d+1 ^ 2 power 
law which is the Fermi liquid behavior well known for 
many decades from weak coupling analysis 0]. It is in- 
structive to realize that the decay of the envelope of the 
overlap function J-(R) for kfR 3> 1 determines the low 
temperature behavior: firstly, only a decaying envelope 
function results in a finite quasiparticlc weig ht Z [hi , 
and secondly, the long-range decay cx Sta,nS/(kfR) 
translates into the power-law behavior of the incoherent 
part of the spectral function and the mobility. 

Here we want to investigate the significance of the sec- 
ond energy scale, the particle-bath interaction U. The 
cross section of the particle is directly related to the 
interaction via the phase shift 5: nkf<7 — M Q T — 
Er sin 2 8 /ird. In phenomenological approaches it was al- 
ways explicitly or implicitly assumed that the interaction 
enters /i(T) only through the cross section. However the 
phase shift may also lead to a subtle interference effect for 
the moving particle. This interference effect is encoded in 
the spatial decay of F(R) (cf. Fig. [j]). It decays anoma- 
lously slow for 8 -> tt/2: F(R) cx 1/ {kjR)^- 1 ^ 2 for 
1 < k f R <C {8t&u8) 2 ^ d ^ ee k f R c whereas for larger R 
the envelope function decays with l/(kfR) d ^ 1 , as stated 
above. For 8 = tt/2, the anomalous decay extends to 
R — > oo. The cross-over length R c corresponds to a 
new low-energy scale E* = En/(Sta,nS) 4 ^ d ~ 1 ^ which ap- 
proaches E* — » with (5 — > 7r/2 and which exists only for 
tan 8 ^> 1, i.e. we consider a strong-coupling anomaly. 
This observation implies that physical quantities will be- 
have differently for T (or w) below or above E* . Specif- 
ically, the power-law behavior above E* is characterized 
by an anomalous exponent — as we will show for /i(T). 
A similar anomalous enhancement of interference effects 
near 8 ~ 7r/2 may be relevant for other fermionic prob- 
lems with dynamical scatterers at distinct spatial posi- 
tions. 
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FIG. 2. ft{T) in the strong coupling regime with 
S — 0.957r/2, uj c ~ €f. The inset displays n(T) on a linear 
scale together with the analytical result of the high-T expan- 
sion. 

We evaluated the mobility in the simplest scheme avail- 
able which still reproduces the exact high T behavior, 
including the non-perturbative logarithmic corrections. 
It also reproduces the Born-approximation results of the 
previous literature which are commonly agreed upon as 
to display the correct power-law behavior for the low- 
est temperature regime. Self-consistency is indispensable 
for a successful evaluation since a diffusing particle ex- 
hibits an entirely different long-time behavior as com- 
pared to a free particle. In order to tackle the mobil- 
ity /ijj(w) = —(l/ioj)(RiRj)(u+io) we have to examine 
a (i?i?)s-correlator which is a path integral evaluated 
with action S (Eq. (0)). For this purpose we introduce a 
general quadratic action 

Sv[A] =S -\ J J^drdr 1 A{t-t') (R{r)-R{r')) 2 (4) 

and expand the correlator to first order in AS = S — Sy: 
(RR}s = (RR) Sv - ((RR)AS) Sv + (RR) Sv (AS} Sv + 
0(AS 2 ). The best possible action Sy[A] to this order re- 
sults from the condition that only the first term survives 
in next-to-leading order, i.e. A(t) has to be found from 
((RR)AS) Sv = (RR)s v ( A S) Sv - And the mobility is 



A(uj + iO,T) 
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where A(lu + iO, T) is the analytical continuation of the 
Fourier transform of A(t). It is straightforward to show 
that this procedure is equivalent to Feynman's variational 
principal (FVP). The early FVP evaluations of the mo- 
bility in the polaron problem did not fully reproduce the 
weak coupling limit. However only few variational pa- 
rameters were used to determine A(r). These inconsis- 
tencies are repaired if the full function A(t) is determined 
self-consistently. Details of this calculation will be pre- 
sented in a more extended article. Fig. ^ displays the 
calculated the static n(T) in d = 3. The phase shift 
is close enough to ir/2 so that the anomalous behavior 
^(T) ss T -3 / 2 extends over two decades. 



The second pronounced strong coupling effect, the gi- 
ant mass renormalization, is illustrated most intuitively 
via a rough estimate based on a Kramers-Kronig rela- 
tion. It relates mass M{u>) and the friction coefficient 
i](lu) as defined by the exact inverse mobility (at T = 0): 
fi(u) + iO)- 1 = r)(w) - iLuM{uj). The friction r)(w) JT|] 
attains its high energy value r]' for frequencies larger 
than the relevant energy scale En (and below Cf). For 
lo <C Eb the phase space for scattering processes ap- 
proaches so that r/(uj) converges algebraically to zero. 
The zero-frequency T = mass enhancement is therefore 
estimated by a step-like rj(oj): AM oc J^f dcur](uj)/uj 2 ~ 

rf/E R = (v'/Er*) ■ {Mo + AM)/M Q . The index o refers 
to unrenormalized quantities. Since the friction coeffi- 
cient is proportional to the cross section of the parti- 
cle, or rather rf/E^ ~ M (kta) 2 where a is the di- 
ameter of the particle, we find from this estimate that 
AM/M - (a/A/) 2 /[l - (a/A/) 2 ] which diverges for 
a ~ Xf. While the divergence is an artifact of this es- 
timate it nevertheless signals that the large number of 
virtual low energy excitations for lu ~ En «6j enhances 
the dynamical mass very effectively 




FIG. 3. Effective mass M(T) for M a = 40m and 
5 — 0.957r/2. The respective number of scattering channels is 
found by comparison with the inset which displays the low-T 
effective mass (the number of channels ranges from 6 to 15). 

The T-dependent dynamical mass has to be calculated 
through the self-consistent scheme outlined above. We 
have to extend it for a large number of scattering channels 
in order to raise a far beyond A 2 . One may include more 
scattering channels by systematically generalizing T(R) 
as has been accomplished by Vladar Q who included 
higher angular momentum channels N up to N = 3. 

However we do not intend to model a specific scatterer. 
In this case, it is more practical to introduce a large num- 
ber of spin channels to investigate this type of strong 
coupling effect in a generic way. The number of scatter- 
ing channels then enters as a multiplicative factor in the 
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exponent of the overlap: Tn{R) = NT(R). Accordingly, 
the inverse high-T mobility scales with a trivial factor TV: 
— N^ 1 — and in leading order of perturbation the- 
ory (J^jv(-R) <C 1) this trivial dependence on N is valid 
down to T — y 0. The self-consistency, however results in 
a tremendous suppression of the low-T mobility (Fig. ||) 
and a giant increase of the zero-frequency effective mass 
(Fig. 0). The numerical analysis shows that 



M 



exp 
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(f) 2 »l (6) 



is an excellent fit (Fig. [|, inset). The factor c depends 
weakly on the phase shift and is approximately c ~ 0.21 
in the perturbative regime and c ~ 0.68 for 5 = 0.95 ir/2. 

An exponential increase of M/M a is indeed expected 
from the requirement for a smooth transition from the 
high-T to the low-T regime: the high-T analysis yields 
for ju(T) a constant plus logarithmic corrections down to 
) @. On the other hand, the low-T 



T, 



Texp{dV/E R 
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behavior is found analytically because self-consistency is 



not required: fJ,(T)- 1 ~ 7i- 2 <5tan<5 • MT (T/Er)^- 1 ^ 2 . 
Assuming that low- and high-T expressions connect 
smoothly at T s , we have to demand that Mt—>o/M > 

2T/T S w cxp^cT/T^) This estimate does not yield 
the same exponent as the full numerical calculation but 
it supports the notion of an exponentially large mass. 



not true for intermediate temperatures: the lower inset 
magnifies this non-universal regime. 

The dramatic rise in the effective mass of large heavy 
particles should be clearly observable for ions in 3 He. 
While the superfluid transition will inhibit the observa- 
tion of the low-temperature effective mass, we expect a 
large increase of the effective mass ex 1/T in the regime 
where a logarithmic increase of the mobility has been 
observed |l(|. The experimental situation is more com- 
plex for muons in metals and will be discussed elsewhere. 
Our predictions are relevant in a regime where the tem- 
perature is low enough so that fi(T) > /i holds where 
fio is the mobility in the absence of a periodic poten- 
tial (which can be determined from fits of the mobility 
to the Kondo-Yamada theory (|J|). Unfortunately, all 
mobility experiments for muons in metals (even in alu- 
minum) seem to be in the opposite regime /i(T) < /xo 
where the effect of the periodic potential determines the 
mobility. Therefore one cither has to reach temperatures 
below lOmK in aluminum or find different metals with 
more mobile muons. 

We acknowledge many helpful discussions with 
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enstein, and especially P. Wolfie. This work was sup- 
ported by the DFG (T.K.). 




10 5 10 4 10 3 10 2 10 1 10° 
(T/(fl c ) * (M(T)/M „) 
FIG. 4. Scaled mobility of the iV-channel model for the 
same parameters as in Fig. 3. The upper inset displays the 
unsealed mobility (the highest mobility corresponds to 6 chan- 
nels, the lowest to 15 channels implying a large cross section). 
The lower inset demonstrates that scaling fails for intermedi- 
ate temperatures. 

Also the mobility depends exponentially on the cross 
section as is displayed in the upper inset of Fig. ^. The 
low-T mobility, however, is explained entirely by the mass 
renormalization. Consequently, all curves of the upper 
inset collapse onto a single line for the lowest T. This is 
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